Nanopattern-stimulated superconductor-insulator transition in thin TiN films 
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We present the results ol the comparative study ol the influence of disorder on transport prop- 
erties in continuous and nanoperforated TiN films. We show that nanopatterning turns a thin 
TiN film into an array of superconducting weak links and stimulates both, the disorder- and mag- 
netic field-driven superconductor-to-insulator transitions, pushing them to lower degree of disorder. 
We find that nanopatterning enhances the role of the two-dimensional Coulomb interaction in the 
system transforming the originally insulating film into a more pronounced insulator. We observe 
magnetoresistance oscillations reflecting collective behaviour of the multiconnected nanopatterned 
superconducting film in the wide range of temperatures and uncover the physical mechanism of 
these oscillations as phase slips in superconducting weak link network. 



That a thin film of the same material can be a su- 
perconductor but can very well turn an insulator, is one 
of the most remarkable aspects of disordered supercon- 
ductors [3-01 ■ The engine driving the transition between 
the superconducting and insulating states is disorder the 
effect of which is two-fold. On the one hand, disorder 
limits the electron diffusion enhancing thus the Coulomb 
electron-electron interaction which competes with the 
Cooper pairing [H [t| . The latter in an interplay with the 
disorder-induced inhomogeneities localizes Cooper pairs 
to form an insulating state, Cooper-pair insulator. A re- 
stricted geometry is critical to effects of disorder - for 
the insulating state to be observed the superconducting 
material is to be thinned down till its thickness d be- 
comes comparable to or smaller than the superconduct- 
ing coherence length £. One of the major experimental 
challenges in these studies remains the optimization of 
material parameters taking it to the closest proximity of 
the direct superconductor-insulator transition and iden- 
tifying the systems that exhibit such a transition at avail- 
able temperatures. In this Letter we meet this challenge 
via creating a metamaterial with the desirable properties, 
the multiconnected thin superconducting film. We show 
that nanopatterning a thin TiN film into a regular sieve- 
like configuration turns it into an array of weak links 
and, therefore, stimulates the direct superconductor-to- 
insulator transition. Depending on the original degree of 
disorder it either suppresses the critical temperature T c , 
or drives the initially superconducting film into an insu- 
lating state, or else, transforms the originally insulating 
film into an even more pronounced insulator. 

As a starting material we have chosen a 5 nm thin 
TiN film which was identical by its parameters to those 
that experienced the superconductor-insulator transition 
after soft plasma etching [1CH13J and which were fully 
characterized by the high resolution electron beam, in- 
frared 14(, and low-temperature scanning tunnelling 



spectroscopy 15j. The smooth, continuous, and uni- 
form TiN film was formed on the Si/Si02 substrate by 
atomic layer deposition. The film had the superconduct- 
ing critical temperature T c = 1.03 K, the diffusion con- 
stant D = 0.3cm 2 /s, and the superconducting coherence 
length £d(0) = 9.3 nm. As a first step, the film was pat- 
terned into the bridges 50 /im wide and 100 /jm long via 
conventional UV lithography. Then, making use of the 
electron lithography and the subsequent plasma etching, 
a square lattice of holes with the diameter ~ 120 nm and 
the period a — 200 nm covering the 50 x 120 /im 2 area, 
was created, see insets in Fig.[TJ The voltage probes VI- 
V2 were designed to fall within the nanopatterned do- 
main of the film, while the probes V3-V4 were placed 
within its continuous (non-perforated) section to measure 
thus the resistance of the original continuous film. The 
perforated section of the film confined between the probes 
V1-V2 contained 50^mxl00/im/(200nm) 2 = 125000 el- 
emental units. To increase the sample sheet resistance, 
it was sequentially treated by an additional soft plasma 
etching two times. We will be referring to the untreated 
sample as to the original one (or as to the state A), dis- 
criminating between the non-perforated (i.e. reference) 
part, rA and the perforated, pA, section of the film. Se- 
quential etching transforms the film into the state B and 
eventually into the state C, and we will be using the pre- 
fixes 'r' and 'p' as in the state A. The room temperature 
resistances are 3.17kf2 for rA, 3.75 kSl for rB, and 4.76 kil 
for rC and 13.87 kfi for pA, 14.45 kQ for pB, and 16.26 kft 
for pC. The ratios of the room temperature resistances of 
the perforated and continuous samples for each respective 
state is about 3-4, reflecting the approximately 3 times 
reduction of the effective cross-section of the sample upon 
perforation. The chosen fabrication procedure allowed 
studying the evolution of electronic transport properties 
in a most controllable way avoiding introduction of addi- 
tional geometric parameters that could vary from sample 
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to sample. The temperature T- and the magnetic field 
-B-dependences of the resistance were measured using the 
standard four-probe dc and low frequency ac techniques. 
The currents were sufficiently small to ensure the linear 
response regime as was verified by direct measurements 
of the current-voltage characteristics I-V. The magnetic 
field was applied perpendicular to the film surface. 

We start discussing our results with the zero magnetic 
field data. Shown in the Fig.[T] are temperature depen- 
dences of the resistance per square (to obtain it the mea- 
sured resistance was divided by two in accordance with 
the sample aspect ratio) for all three states of the perfo- 
rated films [panels (a) and (c)] and for the corresponding 
reference films [panels (b) and (d)]. Both sets demon- 
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FIG. 1. Resistance per square versus temperature plots dis- 
playing superconductor-insulator transitions in the perforated 
films [panels (a) and (c)] and the corresponding reference films 
[panels (b) and (d)]. Panels (a) and (b) employ the linear 
scales of the resistance, whereas in the panels (c) and (d) the 
logarithmic scales are used. The scales of the resistances are 
the same pairwise for the (a)-(b) and (c)-(d) panels. The inset 
to panel (b) sketches the measuring set up (the dimensions 
shown): II and 12 are the current-carrying and the pick up 
leads, respectively, VI and V2 leads are attached within the 
perforated part of the film and V3 and V4 pick up voltage 
from the original non-perforated section of the film. The in- 
set to panel (c) presents scanning electron microscope (SEM) 
image of the part of the TiN film patterned with a square ar- 
ray of holes with the diameter ~ 120 nm and with the center 
spacing of 200 nm. 



strate the superconductor - insulator transitions (SIT). 
To emphasize on the features of the temperature be- 
haviour of the resistance, we present resistance in both, 
linear and logarithmic scales, and the temperature in log- 
arithmic sale since it spans over three orders of magni- 
tude. Resistances of all the films grow upon decreas- 
ing the temperature from the room temperature down, 
with all the superconducting samples showing a non- 
monotonic R(T) behaviour with the pronounced maxi- 
mum preceding the superconducting transition. So, the 
maximal resistance of the as-prepared reference super- 
conducting film rA, is i? max = 5.67kfi. The resistance of 
the film rB achieves i? max = 10.3 M7, and for the nanos- 
tructured film pA the maximal resistance value becomes 
i? max = 39.4kfi at T = 1.64K. On the way down to 
T c = 1.03 K of the rA film the resistance of the pA sample 
shows a noticeable decrease to the value of R — 33.7 kO, 
which is just about of the resistance per square of the 
rA film. Upon further cooling the pA sample exhibits 
the drop in the resistance of about two orders of magni- 
tude. Nevertheless, the pA sample does not transit into 
a global phase coherent superconducting state, remain- 
ing in the resistive state even at lowest temperatures. 
Deferring the discussion of this feature till after the pre- 
sentation of the magnetoresistance data, we stress here 
a striking behaviour of the B-state. While the reference 
sample rB falls to a superconducting state (although it 
occurs at T c — 0.43 K, which is decreased as compared to 
that of rA), the nanostructured part pB appears at the 
insulating state of the SIT. In the state C both, rC and 
pC are insulating. 

We now turn to details of the electronic transport 
properties. Upon cooling down to 10 K all samples ex- 
hibit logarithmic temperature dependence of the conduc- 
tance (see Fig. [2^,), which is well described by the formula 
G(T)/G 00 = Aln(k B Tr/h), where G(T) = l/R a (T) 
is the conductance, Goo = e 2 / {2tt 2 K). This behav- 
ior is in accord with the theory of quantum correc- 
tions for quasi-two-dimensional disordered systems and 
can be attributed to weak localization and repulsive 
electron-electron interaction corrections [l6| . Similar be- 
haviour was observed before in critically disordered TiN 
films [UllZi and in Bi films 0,113. Notably, the high- 
temperature slope is identical for all three samples within 
each sets (with A = 2.6 ± 0.1 for the reference and 
A = 0.85 ± 0.05 for perforated films), irrespectively to 
their low-temperature either superconducting or insulat- 
ing behaviors. The ratio of the factors A in the reference 
and perforated parts of the film is close to 3, as expected, 
due to the geometric reduction of the effective cross sec- 
tions of the conducting channels in nanopatterned struc- 
tures as compared to those in the reference films. This 
evidences that patterning films does not introduce addi- 
tional microscopic disorder. Thus disorder remains the 
same in perforated and continuous films for each state, 
A, B, and C. 
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FIG. 2. (a) Conductance G/G o = 2vr 2 /i/(e 2 ifci) as function 
of temperature. Panel (b): f?n versus 1/T dependences for 
the non-superconducting samples. Both rC and pC show Ar- 
rhenius behaviors, the dashed line depicts the fit to Eq. (f ) 
yielding characteristic temperatures To = 0.63 K and 1.9 K, 
and Ro — 16.6 kfi and 55.6 kfi for samples rC and pC, re- 
spectively. The line bar marks these Tb's emphasizing that 
To's were determined in the domain T < To as proper. At 
lowest temperatures log R(T) turns upwards departing from 
the Arrhenius- for hyperactivated regime [l^. The behav- 
ior of sample pB can not be viewed as the Arrhenius one: 
Linearisation in the interval 1, going up from 0.4 K results 
in To = 0.35 K, and the Arrhenius law cannot be used for 
T > To, while the linearisation in the interval 2 although giv- 
ing acceptable To = 2.21 K, results in an unreasonably low 
Ro — 0.15 kQ. The inset shows a geometry of an elemental 
unit of an array, illustrating the reduction of the cross section 
of the effective conducting channel. 



To analyse the behavior of the non-superconducting 
samples we replot R(T) as function of 1/T in Fig. [2b. 
At low temperatures it is well fitted by the Arrhenius 
formula for both, the reference and the nanopatterned 
films in the state C, evidencing that these samples are 
indeed insulators. The dashed lines correspond to 



R = i? cxp(T /T) 



(1) 



with To = 0.63 K and 1.9 K being the characteristic tem- 
peratures, and i?o = 16.6 kf2 and 55.6 M7 for samples rC 
and pC, respectively. That Ro of the pC is 3.35 times 
larger than that of the reference film, paralleled by the 
ratio of the respective resistances at room temperatures 
that equals to 3.4, gives more support to observation that 
patterning does not change the degree of microscopic dis- 
order. 

Speaking of the characteristic temperatures, one would 
have expected that had To been formed over the micro- 
scopic scale, it would have remained the same in both 



reference and perforated films. Instead, one sees a no- 
ticeable increase in To in a perforated sample. This sug- 
gests that T is built on the macroscopic spatial scales, 
not less then 200 nm and is influenced by changes in geo- 
metric characteristics and the connectivity introduced by 
patterning. Such a behaviour becomes clear, once one re- 
calls that the characteristic energy of an insulating state 
of a two-dimensional Josephson junction array (JJA) is 
&bTo = A c = E c ln(L/6), where E c is the charging energy 
of a single superconducting island, L is the smaller quan- 
tity out of either the electrostatic screening length of JJA 
or its linear dimension, and b is the size of the elemental 
cell of the array 1|1 2(J. Maintaining that the Cooper- 
pair-insulating film comprises a self-induced texture of 
weakly coupled superconducting islands, one would have 
expected here the same behaviour (the size-dependence 
of the characteristic energy was observed in the identi- 
cally prepared TiN films |2l| and InO films (22[). Then 
the upturn in the log R{T) vs. 1/T dependence for the rC 
sample, as well as for the perforated sample pB, occurs 
at the temperature of the charge Berezinskii-Kosterlitz- 
Thouless transition (CBKT) fc B ?CBKT ^ E c ; for the rC 
sample Tcbkt ~ 0.14 K. (In artificially manufactured 
2D JJA the upturn from the behaviour described by the 
Eq. (1) due to CBKT was observed in Refs. [HSl].) This 
yields ln(L/6) = A c /E c = 4.5. In Ref. the quantity 
b w 40 nm was found for critically disordered TiN films 
in the vicinity of the SIT, corresponding to b f=s 4£. We, 
thus, obtain L rs 3.4 /mi, and, since this value is appar- 
ently less then the size of the sample of 50 /im (but at the 
same time well exceeds 200 nm) , L represents the electro- 
static screening length. To crudely estimate the effect of 
perforation we assume the originally square Josephson 
network with the lattice constant f=a 40 nm and notice 
that perforation transforms this network into the new 
one with the 'bond' length equal to 200 nm, each bond 
comprising a series of five junctions. This straightfor- 
wardly gives an estimate for new characteristic energy 
A c = 5A C — 5fJ c ln5 = 2K, which nicely agrees with the 
experimental finding of 1.9 K for the characteristic energy 
for the sample pC. The increase in the characteristic en- 
ergy upon removing the fraction of the junctions of the 
2D JJA was found in 
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Applying magnetic field to nanopatterned films brings 
in competing energy- and spatial scales and thus com- 
mensurability effects which manifest themselves through 
the oscillations in thermodynamics and transport prop- 
erties. Shown in Fig. [3] are the magnetoresistance data 
for the pA sample at relatively weak magnetic field, 
over the range ±0.5 T, which is much less than the up- 
per critical field of the reference film (the reference film 
is close by its parameters to the low resistive sample 
of Ref. [lfj]) B C 2 = 2.8 T. There are eight pronounced 
magnetoresistance (MR) oscillations at low temperatures 
(for each polarity of the field). The main period is 
-Bo = &o/a 2 , where <I>o = irh/e is the superconduct- 
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FIG. 3. The two-dimensional colour map of the resistance 
in the temperature - magnetic field plane for the perforated 
sample pA. The colour scale quantifying magnitudes of the 
resistance is given at the lower right corner. The inset shows 
the experimental data of the magnetoresistance (open cir- 
cles) vs. reduced magnetic field B/Bo, where Bo = &o/a 2 , 
at temperature 0.11 K. In addition to the fundamental dips 
at B n = Bo, where n is an integer, the secondary dips at 
B/Bo = 1/4, 1/3, 2/5, 1/2, 3/5, 2/3, and 3/4 are observed. 
Dash-dotted line corresponds to the periodic dependence of 
the fractional reduction in average Josephson coupling en- 
ergy on the magnetic flux per unit cell derived in the nearest- 
neighbour contour approximation of Ref. 29). The solid line 
passes through the theoretical values calculated for rational 
B/Bo = p/q, with q < 256, on the superconducting net- 
work Eq. (2) accounting for multi-contour phase synchroniza- 
tion and thus exhibiting higher order commensurability ef- 
fects. 



ing flux quantum, corresponding one flux quantum per 
unit cell. Magnification of the resistive curve (see in- 
set in Fig. [3]) reveals a fine structure reflecting collective 
behaviour of the multiconnected superconducting film 
differing it from the behavior of a single superconduct- 
ing loop. One distinguishes additional well-defined dips 
at B/Bq = 1/4, 1/3, 2/5, 1/2, 3/5, 2/3, and 3/4. The 
previous observations reported the B/Bq = 1/2 feature 
in square Josephson junction arrays [251428] . proximity- 
effect junction arrays [29I - I32V perforated films and su- 
perconducting wire networks [331-36] . The "full" set of 



dips was observed in the MR of square J J A [27|, l28[ and 
in proximity- effect junction arrays 32j and in the criti- 
cal temperature variation in square superconducting wire 
network [33j ] . 

A description of the observed modulated MR is based 
on the solution of the Ginzburg-Landau equation on the 
superconducting network in the presence of the magnetic 
field [33[. The wave function ip for the square super- 
conducting network follows the so-called Harper equa- 
tion il: 



Here magnetic field B appears through the flux per pla- 
quette Ba 2 /$o = $/$o = vlli where p and q are relative 
prime numbers. The matrix, corresponding to this equa- 
tion has q x q dimensions. The fields at which resistivity 
exhibits the dips are defined by the boundaries of the 
energy spectrum as function of the magnetic field, i.e. 
one has to find E(§) = max a [e($, a)}. To this end we 
decompose the matrix discriminant into the polinomial 
of e independent on a and the offset function indepen- 
dent on e [38]. The change in the magnetoresistance is 
calculated as AR(B) = Aarccos 2 [i?($)/4], where A is 
a numerical coefficient, which does not depend on the 
magnetic field. The resulting AR(B) behaviour is shown 
by the solid line in the Fig. [3] and demonstrates an excel- 
lent agreement with the data at low temperatures. The 
half-flux dip results mainly from the contribution from 
the currents in the adjacent loops as it was shown in 
Ref. [2{|. As temperature increases the fine structure is 
smeared out, while the oscillations with the main period 
Bo persist till T ~ 0.7 K, where they are easily resolved 
by taking dR/dB. 

The most intriguing aspect of the observed MR oscil- 
lations is an extremely wide temperature region of their 
presence. The MR oscillations in the perforated films 
and/or superconducting wire networks (SWN), were usu- 
ally found, if measured in the linear response regime, in 
the close proximity to T c in the region AT < 0.02T C [33j, 



34 |39It43| ]. More extended temperature re gion s of the 
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ip n+ i + ip n -x + 2 cos(27r$7i + a)ip n = eip r , 



(2) 



MR oscillations were reported in Refs. |3£ 
Juxtaposing the data obtained for various systems and 
inspecting the systems' geometric characteristics we ob- 
serve that as a rule the low-temperature boundary for 
the MR oscillations to appear corresponds to the tem- 
perature at which the ratio w/^(T) < 5, where w 
is the width of the superconducting constriction and 
£d(T) = £ d (0)/y/l-T/T c . This brings to the mind the 
Likharev's result that the weak link cannot accom- 
modate an Abrikosov vortex if w < w c and transforms 
into a Josephson junction, where the critical width w c (L) 
was evaluated as being equal to 4.41£(T) for the long, 
L 3> w, link (L is the link length) near T c . For the 
square link, L — w, w c « 5£ and can become well larger 
for the short weak links. In order to gain the insight into 
the meaning of the boundary separating the dissipative 
behaviour governed by the vortex motion from the resis- 
tive state due to phase slips, let us employ the approach 
used in to evaluate the condition of clustering of vor- 
tex cores. Namely, note that either in the presence of 
the magnetic field or upon passing the current, super- 
conductivity near the edges of the constriction is sup- 
pressed and the Andreev states separated by minigaps 
~ A/ (fc P £) should form. Due to overlap of the wave func- 
tions localized at the opposite edges these states broaden 
as Aexp(— w/£). At the width u> ~ £ln(fc F £), where 
the quasiparticle levels broadening becomes of the or- 
der of the level separation, the constriction turns metal- 
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lie and start to behave as a proximity effect generated 
Josephson junction. The estimate for wq (at the plausi- 
ble values of fc F and £) is in a reasonable agreement with 
the all available experimental data and gives the ratio 
w o/£,(T) — 44-8. Note that the above qualitative consid- 
eration does not use closeness to T c and implies that the 
concept of the critical size of the superconducting weak 
link below which it turns into a Josephson junction can 
be extended to low temperatures. 

The next observation is that in order to exhibit os- 
cillations in the magntoresistance, the SWN should be 
in a resistive state. Thus the second condition deter- 
mining the 'range of observability' of the oscillations is 
that the temperature should be higher than the tem- 
perature of the vortex BKT transition in the SWN, 
Tvbkt- Using the magnitude of the critical current 
observed in our experiment at T — O.IK evaluated 
from the position of the maximum of the dV/dl vs. I 
curve [32], I c = 0.17 /xA, giving I c u\ = 0.68 nA per one 
constriction, one obtains the Josephson coupling energy 
as Ej/k-a = 0.016 K. Therefore, indeed, even the low- 
est temperature of our experiment (0.1K), where the 
MR oscillations are clearly detectable, lie in the inter- 
val T > T VBKT = 7rJ5j/(2fc B ) [Elj]. For this reason the 
Rq(T) dependence for the pA sample remains finite at 
all the experimental temperatures. 

Finally, we briefly discuss the MR at large fields be- 
yond the range of oscillations for the pA sample. At 
fields slightly larger than the critical field B C 2 = 2.8 T of 
the reference film [lfj , we observe a huge peak of the MR 
at some B — B msK , followed by the appreciable negative 
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FIG. 4. The resistance per square versus temperature plots in 
the magnetic fields B = 0, 1.0, 1.5, 1.8, 2.0, 2.5, 3.0 T for the 
perforated sample pA. Left inset: The conductance G/Gqo = 
2iv 2 h/ (e 2 Rn) as a function of the temperature in a logarithmic 
scale in the magnetic fields B — 2.5, 3.0, and 11.0 T. Right 
inset: Magnetoresistance isotherms at temperatures T = 0.1, 
0.15, and 0.4 K. The positions of maxima are B max = 2.9, 3.0, 
and 3.3 T, correspondingly. 



MR, as shown at the right inset in Fig. [4] The value of 
B max shifts towards larger fields upon increase in T. The 
similar shifts of the maxima and the subsequent nega- 
tive MR are observed in the reference film (see Fig. lc of 
Ref. [10|) and are well described by the quantum correc- 
tions to the conductivity due to superconducting fluctua- 
tions [HH . The main panel of Fig.|4]shows the evolution of 
the Rq(T) dependences at different magnetic fields. At 
zero field there is a high temperature, T = 1.64 K, max- 
imum (Fig.QJ. At higher fields, but not too far beyond 
B « 2 T, a second, low temperature maximum appears. 
At yet higher fields, as B approaches and further exceeds 
B C 2 of the reference film, the Rq (T)-dependences become 
monotonically decreasing with the temperature increase. 
Notably, in spite of the fact that the MR exhibits huge 
peak, up to 200 kfi per square at 0.1K, it can be hardly 
viewed as the manifestation of the magnetic-ficld-induced 
superconductor-insulator transition. Indeed, the conduc- 
tance G(T) = l/Rfj(T) at fields around and above the 
peak, fits perfectly to the logarithmic temperature de- 
pendence G(T)/Gqq = A\n(kBTT/h), see left inset in 
Fig. 21 with A monotonically decreasing with the growth 
of the magnetic field, A = 0.49, 0.45, and 0.38 at fields 
2.5, 3.0, and 11.0 T, respectively. This decrease is likely 
to reflect the suppression of the contribution from the 
superconducting fluctuations to conductivity. Moreover, 
we observe that the value of A = 0.38, being multiplied 
by the geometry factor of 3, becomes equal to 1.14, which 
is pretty close A = 1 expected for the contribution to 
the conductivity of disordered metals due to the repul- 
sive electron-electron interaction [l^|. We conclude that 
although this sample does not exactly exhibits the mag- 
netic field driven SIT, it still demonstrates the significant 
amplification of the resistivity growth with the decreas- 
ing temperature (the resistance increases by almost the 
factor of 3.5 in the interval from 1 to 0.1K at B = 3T, 
while in the reference film the corresponding factor is 
about 1.4), and thus can be considered as being put at 
the threshold of the transition from the weak to strong 
localization by nanopatterning. 

In summary, we have shown that nanopatterning disor- 
dered superconducting films pushes the SIT to the lower 
degree of microscopic disorder opening the route to con- 
trol the position of the SIT on the phase diagram. We 
have revealed a wide spectrum of phenomena related to 
periodicity of the phase and the absolute value of the 
superconducting order parameter in a superconducting 
network. 
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